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We show that once the effects of vaUey sphtting and intervaUey scattering are incorporated, 
renormahzation group theory consistently describes the metallic phase in silicon metal-oxide- 
semiconductor field-effect transistors down to the lowest accessible temperatures. 



The two-parameter scaling theory of quantum diffu- 
sion in two dimensions [1 , 2] has been remarkably suc- 
cessful in describing the properties around the metal- 
insulator transition (MIT) in electron systems confined 
to silicon inversion layers (MOSFETs) [3, 4, 5]. The the- 
ory is based on the scaling hypothesis that both the re- 
sistivity and the electron-electron scattering amplitudes 
become scale dependent in a diffusive system due to the 
singular long ranged nature of the diffusive propagators, 
T>{q,uj) = \/{Dq^ + in a disordered medium [6, 7]. 
The predicted scale dependencies calculated using renor- 
mahzation group (RG) theory [1] were recently verified 
experimentally in Ref. [3] without any fitting parameters. 
Since the theory considered the valleys to be degenerate 
and distinct, the experiments were limited to tempera- 
tures larger than the characteristic valley splitting and 
intervalley scattering rate (T > 500 mK). The effects of 
scaling are, however, significant at low temperatures and 
it is therefore important to test the scaling hypothesis at 
much lower temperatures. We show that when the RG 
theory is extended to include valley splitting and inter- 
valley scattering [8] the scaling properties in the metallic 
phase can be described quantitatively down to the lowest 
reliably accessible temperatures, T k, 200 mK. 

The evolution with scale (temperature) of the two- 
parameters, namely, the resistance, p, and the electron- 
electron interaction strength, 72, in the spin-triplet chan- 
nel were discussed in detail for p \ (in units of 
irh/e^) in terms of RG theory in Ref. [1]. (In Fermi- 
liquid notation, 72 is related to the amplitude Fq as 
72 = -Fq°/(1 Fg").) The theory predicts that, while 72 
increases monotonically as the temperature is reduced, 
p behaves non-monotonically, changing from insulating 
behavior {dp/dT < 0) at high temperatures to metal- 
lic behavior [dp/dT > 0) at low temperatures, with the 
crossover occurring when 72 attains the value 72 — 0.45. 
Although the maximum value Pmax occurs at a crossover 
temperature T — Tmax, both of which are sample spe- 
cific and hence non-universal, the two-parameter scal- 



ing theory predicts that the behaviors of p(T) /pmax 
and 72 (T) are universal when plotted as functions of 
^ — Pmax ln(T'max/T)- The abovc predictions, including 
the value of 72 , were verified experimentally in Refs. [1, 3] 
in the temperature range where the two valleys may be 
considered to be degenerate and distinct. 

For n-(OOl) silicon inversion layer the conduction band 
has two almost degenerate valleys located close to the X- 
points in the Brillouin zone. While the sharpness of the 
interface of the inversion layer leads to the splitting, A^,, 
of the two valley bands, the atomic scale irregularities 
found at the interface gives rise to a finite intervalley 
scattering rate, /i/tl [9]. The singularity of the diffusion 
modes, especially those in the valley-triplet sector, are 
cut-off at low frequencies as a result [8, 10]. Hence, the 
specific form of the RG equations, which is sensitive only 
to the number of singular modes, depends on if fc^T is 
greater than or less than the scales or/and h/r^. 

The relevant RG equations for the different tempera- 
ture ranges may be combined as follows [8] : 
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The parameter K accounts for the number of singular 
diffusion modes in each temperature range. For tem- 
peratures T > Ty and T±, where fcs^^u = and 
kBT± ~ fi/T±, the two bands are effectively degener- 
ate and distinct; the constant K in this case is propor- 
tional to the square of the number of valleys, Uy, i.e., 
K = nf, ~ 4 (uy — 2 for silicon). In the temperature 
range T± "^T <Ty, the two bands remain distinct but 
are split and hence each valley contributes independently 
to Ao'(6), i.e., K — Uy — 2. At still lower temperatures 
T < Tl, intervalley scattering mixes the two valleys to 
effectively produce a single valley so that K = 1. 

A few important clarifications regarding the use of 
Eq. (1) are discussed below. First, for the case K = 2, 
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when the bands are split but distinct, it has been shown 
that using a single amplitude 72 to describe the interac- 
tion in all the seven {AK — 1 ) modes is an approximation 
that is valid only if the temperature range T± "^T 
is not too wide [8]. In general, when the bands are split 
certain amplitudes evolve differently from 72, thereby ne- 
cessitating the need to go beyond the two-parameter scal- 
ing description [8, 11]. The deviation is large when the 
RG evolution is allowed to proceed to exponentially large 
scales or T <^ Ty. In our case, however, since Tj_, which 
effectively mixes the two bands, is only a fraction smaller 
than Ty , the deviation of the amplitudes is quickly lim- 
ited by T±. We therefore assume that all the amplitudes 
remain degenerate and contribute equally to p, which 
amounts to taking K = 2 in Eq. (1). 

The second point concerns the weak-localization (WL) 
contribution [12] to Eq. (1). It is seen experimentally that 
the phase breaking rate, h/r^p, saturates at low electron 
densities {n < 10^^ cm^^) for T < 500 K. Correspond- 
ingly, a strong suppression of the WL correction is also 
observed in this regime [13]. These observations are con- 
sistent with our results, as is discussed later. We have 
therefore neglected the weak-localization contribution in 
Eq. (1) when analyzing the cases K = 2 and 1 (these are 
the relevant cases at low temperatures). 

In Ref. [3] it was shown that 72 may be determined ex- 
perimentally by exploiting the 6^ dependence of the mag- 
netoconductance Aa{b) = Aa{B, T) = a{B, T) - cr(0, T) 
in a weak parallel magnetic field b — g^sB /ksT < 1. In 
the weak field limit A(t(6) is given as [14, 15] 

e2 

Aa(5) = -0.091 —if 72 (72 + 1) (2) 

TT/l 

Hence the slope of Ao'(6^) provides a direct measure of 
72, given of course that K is known. 
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FIG. 1; Upper panels: p(T) traces (in units of Tih/e^) 
for three different electron densities, Us — 9.87, 9.58 and 
9.14 X 10^°cm~^. Lower panels: Extracted values of 72(7) 
using Eq. 2 using K — 4, for the same electron densities. (See 
Ref. [3] for further details.) The dashed lines are positioned 
at the critical value 7I = 0.45. Note that the maximum in 
p(T) occurs when 72 attains approximately the value 72. 

In the upper panels in Fig. 1, we plot p(T) at zero mag- 
netic field for three different electron densities. They 



show a characteristic non-monotonic behavior as pre- 
dicted in (1). In the lower panels in Fig. 1 we plot the 
extracted values of 72 using Eq. (2) with K — 4, i.e., 
assuming that the valleys are degenerate and distinct. 
The dashed horizontal line marks the point 72 ~ 0.45 
approximately where p(T) attains its maximum value 
in remarkable agreement with Eq. (1). (At these tem- 
peratures quantum coherence is relevant and its con- 
tribution to weak localization, dp/dS, = n^p^, is to be 
added to Eq. (la).) The results of the comparison be- 
tween theory and experiment are presented in Fig. 2. The 
solid squares (■) are the experimental data points for 
Us = 9.1 X lO^^cm"^, reproduced here from Fig. 1(a). 
The solid lines are the predicted theoretical curves for 
p{T) and 72 (r) with the parameters K — 4, Pmax — 0.4 
and Tmax = 2.3 K. (Here, T^ax is the temperature at 
which p(T) attains its maximum value, p{Tmax) = Pmax-) 
The remarkable agreement between theory and experi- 
ment is especially striking given that the theory has no 
adjustable parameters. 

At temperatures below 0.5 K, the experimentally ex- 
tracted values of 72(7") in Fig. 2(b) seem to saturate with 
further decrease in T. We believe that the saturation 
is an artifact of the analysis related to our assumptions 
that both the valley splitting and the intervalley scatter- 
ing are negligible at the lowest temperatures. As noted 
earlier, the large number of valley modes K = reduces 
to just K = Uy for temperature T± "^T <Ty and to just 
K = 1 ior T <T±. In the following, we recalculate 72 (r) 
taking these considerations into account. 

The experimentally extracted values of 72, using K = 2 
and K = 1, are shown in Fig. 2(b) as diamonds (red ♦) 
and stars (blue -k), respectively. The procedure used to 
extract these values are the same as that used for K — 4, 
namely, by fitting the cr(6^) traces in Fig. 1(b) to Eq. (2) 
using the appropriate K values. We find very favorable 
agreement with theory (solid line) if the crossover scales 
are chosen such that Ty « 0.5 K and Tj, Ki 0.2 K. (Note 
that for these temperatures the WL corrections have not 
been included in Eq. (1) for the reasons discussed ear- 
lier.) These values are in good agreement with earlier 
estimates of Ty [16] and T± [17] obtained at higher den- 
sities employing different methods. We checked by di- 
rect calculation using Eq. (1) that the theoretical values 
of p and 72 are not affected significantly when crossing 
these scales, provided that the WL corrections are not 
included below T < 500K. Deviations from the solution 
for K = 4 taking K = 2 and K = 1 are shown in Fig. 2 
as long (red) and short (blue) dashed lines, respectively. 
As can be seen, the deviations are insignificant (almost 
indiscernible) down to T = 0.2K. 

By comparing with experiments we have extended the 
test of the scaling equations (1) down to the lowest re- 
liably measurable temperatures T w 0.2 K. Concerning 
still lower temperatures, i.e., lower than T — 0.2 K, the 
theory predicts (not shown here) that while p{T) satu- 
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FIG. 2: The result of the comparison between theory (lines) 
and experiment (symbols) for p and 72 are presented in (a) 
and (b), respectively. The parameter K = 4 corresponds to 
the case when the two valleys (n„ = 2) are degenerate, i.e., 
T > T^, where T„ « 0.5K is the estimated valley splitting. 
K — 2 corresponds to the temperature range T < Tv, and 
K = 1 corresponds to the region T < T± « 0.2K where the 
intervalley scattering mixes the two valleys to give one valley. 



rates and then begins to drop again at ultra low temper- 
atures (T < 50 mK), 72(7^) rises fast monotonically for 
K — 1. Further tests of these predictions are in progress. 

To conclude, we have shown that if valley splitting and 
intervalley scattering are incorporated into the RG the- 



ory, the latter quantitatively describes the metallic phase 
down to the lowest readily accessible temperatures. The 
extracted values of intervalley scattering time and val- 
ley splitting are in good agreement with those previously 
obtained at higher densities using different methods. 
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